Numerical Bayesian quantum-state assignment for a three-level quantum system 
II. Average- value data with a constant, a Gaussian-like, and a Slater prior 
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This paper offers examples of concrete numerical applications of Bayesian quantum-state assignment methods 
to a three-level quantum system. The statistical operator assigned on the evidence of various measurement data 
and kinds of prior knowledge is computed partly analytically, partly through numerical integration (in eight 
dimensions) on a computer. The measurement data consist in the average of outcome values of N identical 
von Neumann projective measurements performed on N identically prepared three-level systems. In particular 
the large-Af limit will be considered. Three kinds of prior knowledge are used: one represented by a plausibility 
distribution constant in respect of the convex structure of the set of statistical operators; another one represented 
by a prior studied by Slater, which has been proposed as the natural measure on the set of statistical operators; the 
last prior is represented by a Gaussian-like distribution centred on a pure statistical operator, and thus reflecting 
a situation in which one has useful prior knowledge about the likely preparation of the system. The assigned 
statistical operators obtained with the first two kinds of priors are compared with the one obtained by Jaynes' 
maximum entropy method for the same measurement situation. 

In the companion paper the case of measurement data consisting in absolute frequencies is considered. 



PACS numbers: 03.67.-a,02.50.Cw,02.50.Tt,05.30.-d,02.60.-x 
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1. INTRODUCTION 

In this paper we continue our two-part study [ 1 ] with exam- 
ples of concrete numerical applications of Bayesian quantum- 
state assignment methods to a three-level quantum system. 
Since we will consider the same scenario as in the first pa- 
per, to avoid repeating ourselves we therefore refer the reader 
to the first paper for a more detailed and complete account of 
the motivations, explanations, discussions and references on 
the background, theory, formulas, nomenclature, etc, used in 
this paper. The main difference between the two papers lies 
in the type of measurement data considered. In the first paper 
the measurement data consisted in absolute frequencies of the 
outcomes of N identical von Neumann projective measure- 
ments performed on N identically prepared three-level sys- 
tems. Here we will consider the same measurement situation, 
but the measurement data will instead be in the form of an 
average of values being associated to the measurement out- 
comes, in particular 1, 0, and -1. The statistical operator en- 
coding the average value data and prior knowledge is com- 
puted partly numerically and partly analytically in the limit 
when N — > oo, for a constant, and also for two different kinds 
of a non-constant, prior probability distribution, and different 
average value data. A reason for studying data of this kind, 
other than the obvious one that it may have been given to us 
in this form, is that it constitutes an example of more complex 
data than mere absolute frequencies. It is also interesting to 
study this particular kind of data since it enables us to com- 
pare our assigned statistical operators with those obtained by 
instead using Jaynes' maximum entropy method [2] for the 
same measurement situation. The reason for doing this is that 



we want to investigate whether or not this method could be 
seen as a special case of Bayesian quantum-state assignment, 
and if so, try to find the prior that would lead to the same sta- 
tistical operator as the one obtained by using the maximum 
entropy method. 



2. THE PRESENT STUDY 

In this paper we study data D and prior knowledge / of the 
following kind: 

• The measurement data D consist in the average of 
7Y outcome values of N instances of the same mea- 
surement performed on N identically prepared sys- 
tems. The measurement is represented by the extreme 
positive-operator-valued measure (i.e., non-degenerate 
'von Neumann measurement') having three possi- 
ble distinct outcomes {'1', '2', '3'} represented by 
the eigenprojectors {|1)(1|, |0)(0|, where the 
eigenprojectors are labelled by their associated outcome 
values {1, 0, -1), respectively. We consider the limiting 
case of very large N. 

• Three different kinds of prior knowledge / are used. 
Two of them, 7 co and 7 ga , are the same as those given 



in first paper, i.e. a prior plausibility distribution 
p(p\ Ico) dp = gco(p) dp oc dp, 



(1) 



which is constant in respect of the convex structure of 
the set of statistical operators, in the sense explained 
in HI § 3,4]; and a spherically symmetric, Gaussian- 
like prior distribution 
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P(P\ / g a)dp = g ga (p)dp oc exp 



( 



tr[(p-p) 2 r 



dp, (2) 



2 



centred on the statistical operator p. The latter prior ex- 
presses some kind of knowledge that leads us to assign 
a higher plausibility to regions in the vicinity of p. For 
this prior we consider two examples, when p = |1)(1| 
andp = I0X0I- 1 

The third kind of prior knowledge, 7 s i, is represented by 
the prior plausibility distribution 



p(p\h\)dp = g s i(p)dp 



(detp) M+1 dp, 



(3) 



the so called "Slater prior" for a £/-level system, which 
has been proposed as a candidate for being the appro- 
priate measure on the set of statistical operators. [3] 

The paper is organised as follows: In § [3] we present the 
reasoning leading to the statistical-operator-assignment for- 
mulae in the case of average value data, for finite N and in 
the limit when N —> oo. We arrived at the same formulae (as 
special cases of formulae applicable to generic, not necessar- 
ily quantum-theoretical systems) in a series of papers 0>|5l0]. 
In §|4]we present the particular case studied in this paper and 
give the statistical-operator-assignment formulae in this case, 
introduce the Bloch vector parametrisation, present the cal- 
culations by symmetry arguments and by numerical integra- 
tion, discuss the result and in some cases compare it with that 
obtained by the maximum entropy method. Finally, the last 
section summarises and discusses the main points and results. 



3. STATISTICAL OPERATOR ASSIGNMENT 
3.1. General case 

Again we assume there is a preparation scheme that pro- 
duces quantum systems always in the same 'condition' — the 
same 'state' — where each condition is associated with a sta- 
tistical operator. Suppose we come to know that N measure- 
ments, represented by the N positive-operator-valued mea- 
sures {E* : fx = I,..., rk}, k = l,...,N, are or have been 
performed on N systems for which our knowledge I holds. In 
this paper we will be analysing the case when the data is an 
average of a number of outcome values and it will therefore 
be natural to limit ourselves to the situation when the N mea- 
surements are all instances of the same measurement. Thus, 
for all k = l,...,N, (Ef) = {E^}. 

Let us say that the outcomes ii, ■ ■ ■ , k, ■ ■ ■ Jn are or were 
obtained. Since every outcome is associated to an outcome 
value m,-, the average of all outcome values is 



m = ^ m lk IN. 



(4) 



We will consider the general situation in which the data D 
consists in the knowledge that the average value min N repe- 
titions of the measurement lies in a set T; 



D = [m € T]. 



(5) 



Such kind of data arise when the the measurements is af- 
fected by uncertainties and is moreover "coarse-grained" for 
practical purposes, so that not precise average values are ob- 
tained but rather a region of possible ones. On the evi- 
dence of D we can update the prior plausibility distribution 
g(p)dp := p(p\ I) dp. By the rules of plausibility theory 2 



p{p\DM)dp 



p(D\ p)g(p)dp 
J s p(D\ p)g(p)dp' 



(6) 



where S is the set of all statistical operators. 

The plausibility of obtaining a particular sequence of out- 



comes is 



E/JP) 



1=1 



(7) 



with the convention, here and in the following, that only fac- 
tors with Nj > are to be multiplied over, and where we have 
used that tr{E,p} = p(E,| p) and (Ni,..,N r ) -. N are the ab- 
solute frequencies of appearance of the r possible outcomes 
(naturally, N, > and 2iM = AO. Since the exact order of 
the sequence of outcomes is unimportant and only the abso- 
lute frequencies of appearance N matter, the plausibility of the 
absolute frequencies N inN measurements is 



p(n\ P ) = m f| 



i=i 



[tr{E ; pj ] 

Nil 



(8) 



Define W N as the set of all absolute frequencies Af, for fixed Af 
and r. By the rules of plausibility theory we then have that 



p(D\p)= J] p(D\N A p)p(N\p). 



(9) 



Given that we know Af, we can with certainty tell if Af corre- 
sponds to an average value 



m = 2^ NjTHi/N 



(10) 



that belongs to the set T, and knowledge of the statistical op- 
erator p is here irrelevant. We thus have that p(D\N A p) = 
p(m G T\N) = 1 if iV e (pr and otherwise, where we have 
defined 



(Pr := {A^ e N^| Effm/N e T). 



(11) 



k=\ 



Note that the case p = is equivalent to the case with p = |1)(1|. 



We do not explicitly write the prior knowledge / whenever the statistical 
operator appears on the conditional side of the plausibility; i.e., p(-\ p) := 

P(-\pjy 



3 



Using this together with equations (0 and (O we obtain: 
P (D\ P )=Y jP (N\ P ) = NlY J Y] MEiP]]N ' 

Neipr 



Nil 



(12) 



Netpr i=l 

Inserting this into equation (|6]l we finally obtain: 



p(p\DM)dp 



2 /<n 

iVeiJ>r § 



(13) 



Nil 



■)g(P)dp 



We saw in the first paper that generic knowledge / can be rep- 
resented by or "encoded in" a unique statistical operator: 



J 5 



Pl'= pp(p\l)dp. 



(14) 



The statistical operator encoding the joint knowledge D A / is 
thus given by 



Pdm 



Z /(fl 

Neipr § '=1 



[tr{E iP }] 



(15) 



3.2. Large-W limit 

Let us now summarise some results obtained in [6J] for 
the case of very large N. Consider the general situation 
in which each data set Dn consists in the knowledge that 
the relative frequencies / = (/;) :- (N,/N) lie in a region 
®n — {/I VLifi m i\ e TV}> where is a region in which the 
average values lie (being such that has a non-empty inte- 
rior and its boundary has measure zero in respect of the prior 
plausibility measure). Mathematically we want to see what 
form the state-assignment formulae take in the limit N — > oo. 
Consider a sequence of data sets {Dn}™ =1 with corresponding 
sequences of regions {Yjv}™ =1 and {<Z>iv}^ =I , and assume the 
regions converges (in a topological sense specified in Ja]) to 
regions T ra and (fr^ (the latter also with non-empty interior 
and with boundary of measure zero), respectively. 

Given that the statistical operator is p, the plausibility dis- 
tribution for the outcomes is 

q(p) = (qi(p)) with qi {p) := tr(E ; p). (16) 

In Id] it is shown that 



Further it is also shown that if T M degenerates into a single 
average value in, the expression above becomes 3 

p[p\ in A I] dp oc p(p\ I) 8[Z iqi (p)mi - m] dp. (18) 

This is an intuitively satisfying result, since in the limit when 
iV — > oo we would expect that it is only those statistical opera- 
tors p whose expectation value 2i?!'(P) m i is equal to the mea- 
sured average value m that could have been the case. The data 
single out a set of statistical operators, and these are then given 
weight according to the prior p(p\ I) dp = g(p) dp, specified 
by us. 

By normalising the posterior plausibility distribution in 
equation (fT8l , the assigned statistical operator in equation (PPfl i 
is then given by 



P g(p) S[Ei?/(P) m - m] dp 



Pdm 



I 



(19) 



g(P) 8[£i9/(p) m - m] dp 



4. 



AN EXAMPLE OF STATE ASSIGNMENT FOR A 
THREE-LEVEL SYSTEM 

4.1. Three-level case 



We will now consider the particular case studied in this 
paper. The preparation scheme concerns three-level quan- 
tum systems; the corresponding set of statistical operators 
will be denoted by S3. We are going to consider the case 
when the number of measurements N is very large and in 
the limit goes to infinity. The measurements are all in- 
stances of the same measurement, namely a non-degenerate 
projection-valued measurement (often called 'von Neumann 
measurement'). Thus, for all k = 1, . . ., N, {E® } = {E^} := 
{ 1 1 >< 1 1 , |0><0|, I— 1 ><— 1 1 }, where the projectors, labelled by the 
particular outcome values (m\,m,2,mi) = (1,0,-1) we have 
chosen to consider here, define an orthonormal basis in Hilbert 
space. All relevant operators will, quite naturally and advan- 
tageously, be expressed in this basis. We have for example 
that q^ip) = ti^E^p) = p w , the ^tth diagonal element of p 
in the chosen basis. As data we are given that the average of 
the measurement outcome values is in (more precisely in the 
sense that i> ro degenerates into a single average value in). 



4.2. Bloch vector parametrisation and symmetries 

We will be using the same parametrisation of the statistical 
operators as in the companion paper, i.e. in terms of Bloch 



p(p\D N A I) dp cc 



|0, iff(p)e (Doo, 

[p(p\I)dp, if q(p)eO m , 

as N 



(17) 



3 Note that we have here, with abuse of notation, written p[p\ m A /] instead 
of the more correct form p[p\ (fh = m*) A /], to avoid introducing another 
variable in* for the average value data. 
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vectors x. For a three-level system the Bloch vector expansion 
of a statistical operator p(x) is given by: 



1 „ 1 



PW = 3/3 + , v/A ; . 



where 



Xi = tr{A; p(x)} = (Ai)p( X )- 



(20) 



(21) 



The Gell-Mann operators A , are Hermitian and can therefore 
be regarded as observables. Note that our von Neumann mea- 
surement corresponds to the observable 



a 3 s|i>UI + oioxoi -i-ix-ii. 



(22) 



Hence, given a statistical operator p(x), the following holds 
for the expectation value of the outcome values {1, 0, -1} for 
this particular measurement: 

<As> P (jt) = T Ji qi{x)m i = p n (x) - P33W = x 3 . (23) 

Equation (TT~8T > thus becomes 

p[x\ in A /] dx oc g(*) 6(^3 - m) dx, (24) 

and the assigned statistical operator in equation dT9b assumes 
the form 



p(x) g(x) 8(x3 - m) dx 



PmM 



(25) 



(jc) 6(x3 - m) dx 



where Bg is the set of all three-level Bloch vectors. This can be 
rewritten in a form especially suited for numerical integration 
by computer, which we shall use hereafter: 



p(x)g(x)8(x3 - m) Xb s (x) dx 



PmM 



J g(*)5(x3 - m) xb s (x) dx 

c 8 



(26) 



where Xr s (x) is the characteristic function of the set Bs c 
Cs := [-1, l] 7 x [ — -nj> -^]- Using the Bloch vector expansion 
in equation (|20b we see that by computing the following set of 
integrals we have determined p mAI : 

U\m, I] := J Xi g(x) S(x 3 - in) Xb s (x) dx, (27) 

where i e {1, .., 8), and 

Z[m, I] := J g(x) S(x 3 - in) Xb s (x) dx, (28) 



where the dependence of the average value and prior knowl- 
edge is indicated within brackets. The assigned statistical op- 
erator will then given by 



3 24-f 

1=1 



1 Li[m, I] 

2 Zj Z[ih,I] ' 



(29) 



One sees directly from equations d27l i and ( |28] | that 
L3 [in, I] /Z[m, I] — in (Z[m, I] can never vanish, its integrand 
being positive and never identically naught). 

For the same reasons already accounted for in the first paper 
we will not try to determine p jSA/ exactly, but also here com- 
pute it with a combination of symmetry considerations of Bs 
and numerical integration. For all three kinds of prior knowl- 
edge considered in this paper the same symmetry arguments 
used in the companion paper also holds here, so again we have 
that Li[m, I]/Z[fh, I] = for all i # 3, 8 and any average value 
-1 < m < 1. The assigned Bloch vector is thus given by 
(0, 0, in, 0, 0, 0, 0, Lg[w, Z]/Z[m, /]). This means that p AA/ lies 
in the (X3, jcg)-plane and it has, in the chosen eigenbasis, the 
diagonal matrix form 



PmM ~ 



2 V3Z[m,/] 





V3Z[m,/] 








m , Lg[m,/] 
2 2 V3Z[m,/] 



(30) 



4.3. Numerical integration, results and the maximum entropy 
method 

We have used numerical integration 4 to compute 
L&[ih, I] /Z[fh, I] for different prior knowledge and dif- 
ferent values of iit. The result for a constant prior density 
is shown in figure Q] where the blue curve (with bars indi- 
cating the numerical-integration uncertainties) is the Bloch 
vector corresponding to p„-, A / plotted for different values of 
JC3 = in. 5 

It is interesting to compare p, ?lA/ with the statistical opera- 
tor obtained by the maximum entropy method [2] for the mea- 
surement situation we are considering here. Given the expec- 
tation value (M) of a Hermitian operator M, corresponding to 
an observable M, the maximum entropy method assigns the 



4 Using quasi Monte Carlo-integration in Mathematica 5.2 on a PC (Pen- 
tium 4 processor, 3 GHz). The computation times are given in figures [T] 
to [4] and for more details on the numerical integration we again refer the 
reader to the companion paper Q||. 

5 Note that we have for all three kinds of priors considered in this pa- 
per computed Lg[m, I]/Z[m, I] only for non-negative values of X3 = fit, 
since by using the symmetry operation (xi,X2,X3,X4,Xs,x^,Xj,Xg) h» 
(X6,xj,-X3,X4,-X5,xi,X2,xs) one can show that Lg[m, r]/Z[m, I] is in- 
variant under a sign change of m. Further, we have not computed 
Lt,[m, r]/Z[m, /] for fh = ±1, since it follows from (J eq. 17] that 
Lg[m, I]/Z[m, /] = 1/ V3 is the only possibility in this case (which one 
also realises by looking at the figures). 
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statistical operator to the system that maximises the von Neu- 
mann entropy S := -trjplnp] and satisfies the constraint 
tr{pM) = (M). Having obtained an average value M from 
many instances of the same measurement performed on iden- 
tically prepared systems, one conventionally sets (M) = M. 

In our case the operator M would be identified as the Her- 
mitian operator A3 and M as m. Hence the maximum entropy 
method corresponds here to an assignment of the statistical 
operator that maximises S among all statistical operators sat- 
isfying (A3) = m, and this statistical operator is given by 



-fi(m) A 3 



Pme '■- 



tr{e 



-H(m)A 3 i 



(31) 



where 



\i(m) := In 



-m + V4 - 3 m 2 
2(m+ 1) 



(32) 



This could be compared with the statistical operator p, ?lA/ ob- 
tained by instead using Bayesian quantum-state assignment, 
and expressed in general form as in equation (|25] l it is seen to 
instead be given by a weighted sum, with weight g(x) dx, of 
all statistical operators with (A3) = = in. 

In the case of a constant prior one sees from figure Q] that 
PmM is in general different from p ME (the red curve [without 
bars]). This means for instance that, if the maximum entropy 
method is a special case of Bayesian quantum-state assign- 
ment, the statistical operator obtained by the former method 
corresponds to a non-constant prior probability distribution 
g(x) dx on Bg in the latter method. This conclusion in itself 
is perhaps not so surprising, but it raises an interesting ques- 
tion: Does there exist a (non-constant on B%) prior distribution 
g(p)dp that one with Bayesian quantum-state assignment in 
general obtains the same assigned statistical operator as with 
the maximum entropy method? 

A strong candidate is the "Bures prior" which has been pro- 
posed as the natural measure on the set of all statistical opera- 
tors (see e.g. J3>d,H 13 HI]), but unfortunately it turns out to 
be difficult to do numerical integrations on it due to its compli- 
cated functional form, so we have not computed the assigned 
statistical operator in this case. Another interesting candidate 
is the "Slater prior" [3J], which have also been suggested to 
be the natural measure on the set of all statistical operators, 
and the computed assigned statistical operator in this case is 
shown in figure [2] One can see directly from the figure that 
although it is similar to the curve obtained by the maximum 
entropy method, we have found them to differ. 
The computed assigned statistical operators for the Gaussian- 
like prior, centred on the projectors p = |1)(1| and p = |0}(0| 
with "breadth" s = 1/4, are shown in figures [3] and [4] respec- 
tively. Apart from being symmetric under a sign change of m, 
as already have been noted in footnote [5] one can also show 
that L%\m, I g3 \/Z[m, 7 ga ] does not depend on the xj -coordinate 
of the statistical operator the prior is centred on. 



5. CONCLUSIONS 

This was the second paper in a two-part study where the 
Bayesian quantum-state assignment methods has been applied 
to a three-level system, showing that the numerical implemen- 
tation is possible and simple in principle. This paper should 
not only be of theoretical interest but also be of use to ex- 
perimentalists involved in state estimation. We have anal- 
ysed the situation where we are given the average of outcome 
values from N repetitions of identical von Neumann projec- 
tive measurements performed on N identically prepared three- 
level systems, when the number of repetitions N becomes very 
large. From this measurement data together with different 
kinds of prior knowledge of the preparation, a statistical oper- 
ator can be assigned to the system. By a combination of sym- 
metry arguments and numerical integration we computed the 
assigned statistical operator for different average values and 
for a constant, and also for two examples of a non-constant, 
prior probability distribution. 

The results were also compared with that obtained by the max- 
imum entropy method. An interesting question is whether 
there exists a prior probability distribution that gives rise to 
an assigned statistical operator which is in general identical 
to the one given by the maximum entropy method, i.e. if the 
maximum entropy method could be seen as a special case of 
Bayesian quantum-state assignment? In the case of a constant 
and a "Slater prior" on the Bloch vector space of a three-level 
system we saw that the assigned statistical operator did not 
agree with the one given by the maximum entropy method. It 
would therefore be interesting to try other kinds of priors, in 
particular "special" priors like the Bures one. 
The generalisation of the present study to data involving dif- 
ferent kinds of measurement is straightforward. Of course, in 
the general case one has to numerically determine a greater 
number of parameters (the Lj[m, I]) and therefore compute a 
greater number of integrals. 

Post scriptum: During the preparation of this manuscript, 
P. Slater kindly informed us that some of the integrals numer- 
ically computed here and in the previous paper can in fact be 
calculated analytically, using cylindrical algebraic decompo- 
sition |[l2 , 13 , 14 , 15 , 3 with a parametrisation introduced by 



Bloore~ jfl7ll : cf. Slater Ill8ll . This is true, e.g., for the integrals 
involving the constant and Slater's priors. By this method 
Slater has also proven the exact validity of eq. (52) of our pre- 
vious paper [1]. We plan to use and discuss more extensively 
this method in later versions of these papers. 



Acknowledgements 

We cordially thank P. Slater for introducing us to cylindrical 
algebraic decomposition and showing how it can be applied to 
the integrals considered in our papers. AM thanks Professor 
Anders Karlsson for encouragement. PM thanks Louise for 
continuous and invaluable support, and the staff of the KTH 



6 



Biblioteket for their irreplaceable work. 

(Note: 'arxiv eprints' are located at http : //arxiv . org/ ) 



[1] A. Mansson, P. G. L. Porta Mana, and G. Bjork, Numeri- 
cal Bayesian quantum-state assignment for a three-level quan- 
tum system. I. Absolute-frequency data with a constant and a 
Gaussian-like prior (2006), arxiv eprint quant-ph/0612105 

[2] E. T. Jaynes, Information theory and statistical me- 
chanics. II, Phys. Rev. 108(2), 171- 190 (1957), 
http://bayes.wustl.edu/etj/nodel.html see also 
Information theory and statistical mechan- 
ics, Phys. Rev. 106(4), 62 0-630 (1957), 
http://bayes.wustl . edu/et j /node 1 .html 

[3] P. B. Slater, Reformulation for arbitrary mixed states of Jones' 
Bayes estimation of pure states, Physica A 214(4), 584-604 
(1995). 

[4] P. G. L. Porta Mana, A. Mansson, and G. Bjork, From "plausi- 
bilities of plausibilities " to state-assignment methods: I. "Plau- 
sibilities of plausibilities " : an approach through circumstances 
(2006), arxiv eprint quant -ph/0607111 

[5] P. G. L. Porta Mana, A. Mansson, and G. Bjork, From "plausi- 
bilities of plausibilities " to state-assignment methods: II. In- 
duction and a challenge to de Finetti's theorem (2006), in 
preparation. 

[6] P. G. L. Porta Mana, A. Mansson, and G. Bjork, From "plau- 
sibilities of plausibilities" to state-assignment methods: III. In- 
terpretation of "state " and state-assignment methods (2006), in 
preparation. 

[7] M. S. Byrd and P. B. Slater, Bures measures over the spaces 
of two- and three-dimensional density matrices, Phys. Lett. A 
283(3^), 152-156 (2001), arxiv eprint quant-ph/0004©55 



[8] P. B. Slater, Bures geometry of the three-level quantum sys- 
tems, J. Geom. Phys. 39(3), 207-216 (2001), arxiv eprint 
|quant-p h/0008069 , see also d]. 
[9] P. B. Slater, Bures geometry of the three-level quantum systems. 
II (2001), arxiv eprint math-ph/0102032 ; see also OJ. 

[10] P. B. Slater, Hall normalization constants for the Bures volumes 
of the n-state quantum systems, J. Phys. A 32(47), 8231-8246 
(1999), arxiv eprint quant-ph/9904101 

[11] P. B. Slater, Applications of quantum and classical Fisher infor- 
mation to two-level complex and quaternionic and three-level 
complex systems, J. Math. Phys. 37(6), 2682-2693 (1996). 

[12] D. S. Arnon, G. E. Collins, and S. McCallum, Cylindrical alge- 
braic decomposition I: The basic algorithm, SIAM J. Comput. 
13(4), 865-877 (1984), see also GJ]. 

[13] J. H. Davenport, Y. Siret, and E. Tournier, Computer Algebra: 
Systems and Algorithms for Algebraic Computation (Academic 
Press, London, 1987/1993), 2nd ed., transl. by A. Davenport 
and J. H. Davenport; first publ. in French 1987. 

[14] B. Mishra, Algorithmic Algebra (Springer- Verlag, New York, 
1993). 

[15] M. Jirstrand, Cylindrical algebraic decomposition 
— an introduction, Tech. Rep. LiTH-ISY-R-1807, 
Linkoping University, Linkoping, Sweden (1995), 
http : //www . control . isy . liu . se/publications/ doc?id=164 

[16] C. W. Brown, Simple CAD construction and its applications, J. 
Symbolic Computation 31(5), 521-547 (2001). 

[17] F. J. Bloore, Geometrical description of the convex sets of states 
for systems with spin-^ andspin-l, J. Phys. A 9(12), 2059-2067 
(1976). 

[18] P. B. Slater, Two-qubit separability probabilities and beta func- 
tions (2006), arxiv eprint quant-ph/0609006 

[19] D. S. Arnon, G. E. Collins, and S. McCallum, Cylindrical alge- 
braic decomposition II: An adjacency algorithm for the plane, 
SIAM J. Comput. 13(4), 878-889 (1984), see also El- 



7 




Figure 1: Bloch vectors of the assigned statistical operator for prior knowledge / co , computed by numerical integration for different average 
values in = X3 = (A3) (connected by the blue curve [with bars]). The red curve (without bars) is the statistical operator given by the maximum 
entropy method also as a function of the average value in = x 3 = (A 3 ). The large triangle is the two-dimensional cross section of the set 
B 8 along the plane Ox 3 x$. The maximum numerical-integration uncertainty in the x% component is ±0.01. Note that only the ten points for 
< in < 0.9 have been determined by numerical integration, since the nine points for -0.9 < in < -0.1 can be exactly determined from 
the former by symmetry arguments. The endpoints corresponding torn = ±1 were set manually, since x g = 1/ V3 is the only possibility in 
this case. Within the given uncertainties, numerical computations yielded the exact results. The computation was done on a PC (Pentium 4 
processor, 3 GHz) and the computation time was 15 min. 
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Figure 2: Bloch vectors of the assigned statistical operator for prior knowledge / s i, computed by numerical integration for different average 
values m = x 3 = (/\ 3 ) (connected by the blue curve [with bars]). The red curve (without bars) is the statistical operator given by the maximum 
entropy method also as a function of the average value in = x 3 = (A 3 ). The large triangle is the two-dimensional cross section of the set 
B 8 along the plane 0x 3 x g . The maximum numerical-integration uncertainty in the x 8 component is ±0.02. Note that only the ten points for 
< m < 0.9 have been determined by numerical integration, since the nine points for -0.9 < in < -0.1 can be exactly determined from 
the former by symmetry arguments. The endpoints corresponding torn = ±1 were set manually, since x g = 1/ V3 is the only possibility in 
this case. Within the given uncertainties, numerical computations yielded the exact results. The computation was done on a PC (Pentium 4 
processor, 3 GHz) and the computation time was 250 min. 
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Figure 3: Bloch vectors of the assigned statistical operator for prior knowledge / ga , computed by numerical integration for different average 
values m = X3 = (A3) (connected by the curve). The large triangle is the two-dimensional cross section of the set B 8 along the plane Ox^Xg. 
The prior knowledge is represented by a Gaussian-like distribution of "breadth" s = 1/4 centred on the pure statistical operator |1)(1|. The 
small circular arc is the locus of the Bloch vectors (on the plane) at a distance \x — Jt\ = s from the vector x := (0,0, 1,0,0,0,0, 1/ V3) 
corresponding to the statistical operator The numerical-integration uncertainty in the xg component is ±0.016. Note that only the ten 

points for < in < 0.9 have been determined by numerical integration, since the nine points for -0.9 < m < -0. 1 can be exactly determined 
from the former by symmetry arguments. The endpoints corresponding to m = ±1 were set manually, since Xg = 1/ V3 is the only possibility 
in this case. Within the given uncertainties, numerical computations yielded the exact results. The computation was done on a PC (Pentium 4 
processor, 3 GHz) and the computation time was 30 min. 
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Figure 4: Bloch vectors of the assigned statistical operator for prior knowledge / ga , computed by numerical integration for different average 
values m = X3 = (A3) (connected by the curve). The large triangle is the two-dimensional cross section of the set B 8 along the plane Ox^Xg. 
The prior knowledge is represented by a Gaussian-like distribution of "breadth" s = 1/4 centred on the pure statistical operator |0)(0|. The 
small circular arc is the locus of the Bloch vectors (on the plane) at a distance \x — Jt\ = s from the vector x := (0, 0, 0, 0, 0, 0, 0, -2/ V3) 
corresponding to the statistical operator |0)(0|. The numerical-integration uncertainty in the xg component is ±0.02. Note that only the ten 
points for < in < 0.9 have been determined by numerical integration, since the nine points for -0.9 < m < -0. 1 can be exactly determined 
from the former by symmetry arguments. The endpoints corresponding to m = ±1 were set manually, since Xg = 1/ V3 is the only possibility 
in this case. Within the given uncertainties, numerical computations yielded the exact results. The computation was done on a PC (Pentium 4 
processor, 3 GHz) and the computation time was 425 min. 



